Revival times at quantum phase transitions 
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The concept of quantum revivals is extended to many-body systems and the implications of 
traversing a quantum phase transition are explored. By analyzing two different models, the vibron 
model for the bending of polyatomic molecules and the Dicke model for a quantum radiation field 
interacting with a system of two-level atoms, we show evidence of revival behavior for wave packets 
centered around energy levels as low as the fundamental state. Away from criticality, revival times 
exhibit smooth, nonsingular behavior, and are proportional to the system size. Upon approaching 
a quantum critical point, they diverge as a power law and scale with the system size, although the 
scaling is no longer linear. 



PACS numbers: 05.70.Jk,42.50.Md 
I. INTRODUCTION 



The concept of a phase transition can be extended to 
zero absolute temperatures, when thermal fluctuations 
cease. Quantum fluctuations then take over and the sys- 
tem may undergo a quantum phase transition (QPT), 
which reflects in a dramatic change in its physical prop- 
erties as illustrated by the dependence of many several 
observables on a suitable control parameter (other than 
the temperature, say A) that determines the amplitude 
of quantum fluctuations [JJ. It is a generic feature that 
at second order QPTs the system will exhibit diverging 
quantities as the transition is approached. An important 
one is the correlation length, which, following notions and 
a nomenclature borrowed from classical critical phenom- 
ena, diverges as £ ~ |A — X c \~ v - Similarly, other quanti- 
ties such as specific heats or susceptibilities diverge in the 
same manner, defining new critical exponents which are 
related through scaling relations. This is all well known 
and to date examples abound of quantum critical fea- 
tures of classical observables brought about by quantum 
fluctuations alone. Much less studied, however, is the 
influence of QPTs on purely quantal properties, lacking 
a classical counterpart. Here, we address the problem of 
wave-packet revivals in systems exhibiting QPTs. 

The long-time evolution of propagating quantum wave- 
packets may lead to unexpected periodic behavior. Ini- 
tially, the wave packets evolve quasiclassically and oscil- 
late with a classical period Xci, but eventually spread 
out and collapse. At later times, multiples of the revival 
time Tr, wave packets regain their initial form and be- 
have quasiclassically again. The classical period and the 
revival time of wave packet evolution are embodied in 
the first coefficients of the Taylor expansion of the en- 
ergy spectrum Ek around the energy Ek corresponding 



to the peak of the initial wave packet, 

E k = E ko +E' ko {k-k Q ) + ^{k-k ) 2 + ^{k-k f + - ■ ■ , 

(1) 

the first-, second- and third-order terms in the expansion 
providing the classical period of motion Tq\ = 2%/\E' k |, 
the quantum revival scale time Tr = ATi/\E' k ' |, and the 
so-called super-revival time Tsr = 12^/1^^," |, respec- 
tively. Revivals have received considerable attention over 
the last decades [2J. Both experimental and theoretical 
progress was made in, among others, Rydberg atoms, 
molecular vibrational states, electric currents in graphene 
or Bose-Einstein condensates [3HS]- Recently, methods 
for isotope separation [TIT], number factorization [TJJ as 
well as for wave-packet control [T^l - TH] have been put 
forward that are based on revival phenomena, and the 
presence of effective multi-body interactions in a system 
of ultracold bosonic atoms in a three-dimensional optical 
lattice was identified in time-resolved traces of quantum 
phase revivals [15] . Interestingly, the collapse and revival 
dynamics of ultracold atoms in optical lattices have been 
investigated and shown to be strongly sensitive to the 
initial many-body ground state |16j . 

In this article, we extend the concept of quantum re- 
vivals to many-body systems and explore the implica- 
tions of traversing a quantum phase transition. By an- 
alyzing two different models, the U(3) vibron model for 
the bending dynamics of molecules and the Dicke model 
of two-level atoms interacting with a one-mode radia- 
tion field, we show that, as a consequence of the squeez- 
ing of the energy levels around the ground-state energy 
at the critical point, revival times of wave packets cen- 
tered around quantum numbers as low as the fundamen- 
tal state diverge as a power law upon approaching a quan- 
tum critical point. Interestingly, we find that the revival 
times diverge sufficiently close to the critical point even 
at finite system-sizes. 
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II. THE U(3) VIBRON MODEL 



The U(3) vibron model has been successfully applied 
to study the bending dynamics of linear polyatomic 
molecules [T7]. Its Hamiltonian is constructed as a com- 
bination of invariant operators associated with the sub- 
algebras of U(3) (see [17] for further details), and it reads 



H = (l- X )k 



N-V 



-P, 



(2) 



where < x < 1 and P is the so called pairing opera- 
tor (see below). Here N and k are the total number of 
bound states and the number of vibrational quanta. The 
associated base has elements \k, l), I being the value of 
the angular momentum along the z axis perpendicular to 
the plane of vibrations. Their allowed values for a given 
N axe k = 0,1,2, ... ,N and I = ±k, ±(k - 2), . . . , ±1 
or for k odd or even. In what follows we shall restrict 
ourselves to zero vibrational angular momentum, I = 0. 

Four cases can be distinguished in terms of the value 
of the control parameter x- bending vibrations of rigidly 
linear molecules (x = 0); the quasilinear case (0 < x < 
0.2), (including a quantum critical point at x = 0-2; see 
below); the quasibent case (0.2 < x < l)j and the rigidly 
bent case (x = !)■ The Hamiltonian is U(2) invariant for 
X = and SO(3) invariant for x = 1, and in these cases 
analytical solutions for the spectra exist [T8] : 



E k { X = 0) = 2k, 
E k (.X = 1) 



N — 1 



N 



k-k 2 



(3) 
(4) 



The classical and revival times can be computed from 
the expressions above to give the exact, limiting values 
Tc\(x = 0) = 7r and Tr(0) = oo on the one hand, and 
Tab = 1) = (N- 1)tt/(2JV + 1 - 4fc ) and T R (1) = 
(N — l)7r/2 on the other. The divergence of Tr is in 
accordance with the harmonic nature of the Hamiltonian 
when x = 0, while Tr(x — 1) being linear in N provides 
an example of the scaling of revival times in a many-body 
system, a scaling feature that is absent in the classical 
period. 

For general < x < 1 the spectra are no longer an- 
alytically accessible and one has to resort to numerical 
methods. To diagonalize the Hamiltonian Eq. |2]) we use 
the matrix elements of P = N(N + 1) - W 2 [IS], 



(k 2 ,l\W 2 \fa,l) = a5 kl ,k 2 + bSk!-2,k 2 + c5kt+2,k 2 i ( 5 ) 

(N - fa)(fa + 2) + (N-fa + l)fa + I 2 -, 
-^{N -k!+ 2)(N -h + l)(fc! + l)(k ± - 0, 
c = -y/(N - fci)(JV - fa - l)(fci + I + 2)(fci -1 + 2). 

(6) 

We find the general feature (already reported in [T§] ) 
that on increasing the system size TV, k/N becomes a 



with 

a - 
b 




FIG. 1. (Color online) Time dependence of the modulus of 
the autocorrelation function |A(t)| for wave packets initially 
centered around k /N = 0.25 with N = 1000 (bottom), 2000 
(center), and 4000 (top). Other parameter values are \ = 0.5 
and cr = 2. Tr ~ 1024 denotes the revival time for N = 
1000, and Tqi — 192 the classical period. Times are given in 
dimensionless units. 



quasi continuous variable x and that the combination 
Ek ^ N converges to a common spectrum e x . This property 
alone implies that the revival times for a given x scale 



with N because 



d 2 e x 
dx 2 



N- 



. d 2 E k 
dk 2 ' 



(7) 



For the purpose of illustrating this scaling, we construct 
initial wave packets as the linear combination 



|*(t = o)> = 53c fc |M = o) 



(8) 



with Gaussian coefficients c k oc exp[— (k — ko) 2 /a], a = 2 
and centered around Xq = k^/N = 0.25. Centering the 
wave-packets around a common xq guarantees that the 
average energies are common in turn |20) . At later times, 
of course, 



|M/(t))=^c fc |M = 0}e 



-iE k t 



(9) 



Figure [T] shows the time evolution of the modulus of the 
autocorrelation function, A(t) = (^f (Q)\^! (t)) , which is 
the overlap between the initial and the time-evolving 
wave packet, of three wave-packets corresponding to 
X = 0.5, N = 4000, 2000, and 1000. Given an initial 
state, \A(t)\ decreases in time and the occurrence of re- 
vivals is reflected in its returning to its initial value of 
unity. 

It can be clearly appreciated from the figure that 
T R (N = 4000) ~ 2T R (N = 2000) ~ 4T R (N = 1000). 
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FIG. 2. (Color online) Time dependence of |j4(t)| for a wave- 
packet centered around the ground state and system param- 
eters x = 0.5, TV = 2000 and a = 2. T R ~ 2853.5 denotes the 
revival time. Time is given in dimensionless units. 



FIG. 3. (Color online) Revival and super-revival (inset) times 
at criticality as a function of \ f° r wave-packets centered 
around the ground state and TV = 1000. The red circles 
and the blue squares correspond to, respectively, x < Xc and 
X > Xc- Times are given in dimensionless units. 



For consistency, wc have also verified that these revival 
times match those obtained by evaluating 4n/\E'£ |, the 
second derivatives being simply computed through the 
numerical approximation E^+i + Ek-i — 2Ek- The figure 
also shows the classical period Tq\ ~ 192, which does not 
scale with TV. 

Revivals are also observed at energy levels as low as 
the ground state. As an example, Fig. [2] shows the time 
development of a wave packet with Gaussian-distributed 
population for x — 0.5, centered around the ground 
state and with a — 2. The estimated revival time from 
the figure is approximately 2850, to be compared with 
T R /2 = 27r/|££ o=0 | = 2853.5. In this case, E% is com- 
puted by fitting the first three levels of the spectrum to 
a parabola and then taking k = 0. The results thus 
obtained are in perfect agreement with those observed 
by monitoring the wave-packets time evolution. Other 
values of x, both above and below 0.2, yield again the 
behavior Tr ~ TV. 

Next, we study how the revival times are affected by 
the presence of the quantum critical point. At x = 
Xc = 0.2 this system undergoes a second-order quantum 
phase transition in the thermodynamic limit, TV — > oo, 
between two phases displaying anharmonicitics of oppo- 
site signs. To compute the classical period and the re- 
vival time, notice that the former is simply related to the 
energy gap, A = E\ — E , by Tq\ ~ A -1 . Since this 
can be evaluated in the thermodynamic limit through 
A = ^{5 X ~l){l + 3x) (valid for X > Xc H5), this im- 
plies that the classical period of wave-packets centered 
around the ground state diverges as the critical point 
is approached as Tc\ ~ (x ~ Xc) 1 ^ 2 - For the sake of 
consistency, we have verified this prediction numerically. 
As regards Tr, a divergent behavior is also expected due 



to the squeezing of the energy levels around the ground 
state energy at the critical point, what leads to the time 
evolution of the wave-packet being basically controlled 
by a single autostate. In fact, we find Tr ~ (x~Xc) 1 at 
fixed TV = 1000 and fc = 0, with X c = 0.205 907075(3) 
(see Fig. |3j. Supporting this image is the divergence 
of the super-revival time, T S r, as shown in the inset of 
Fig. [3] with, in this case, X c = 0.203 9044(2) and the 
same scaling exponent -1. Notice, that all time scales 
diverge at close but different x> on ly to coincide at \c 
in the thermodynamic limit [for TV = 1000, Tq\ peaks at 
Xc = 0.205 305(5) rather than diverging]. 

Another interesting aspect is the dependence of Tr and 
Xci on the system-size TV for a given k. Previous results 
for the energy gap yield the scaling form A ~ TV -1 / 3 (see 
[PE] and references therein) , which is in perfect agreement 
with our own, T C i(A:o = 0) - TV 1 / 3 . Moreover, we find 
this same scaling behavior for wave-packets peaked at ar- 
bitrarily high k, i.e. T cl (fc > 0) - TV 1 / 3 . Turning to T R , 
we find that for any system-size TV there always exists a 
X, the closer to X c the larger TV, such that T^ (k = 0) 
diverges. This is illustrated in the inset of Fig. [4] which 
shows the behavior of Tr for TV = 1000 in the vicinity of 
Xc- Similar divergences are found generically for any TV, 
including values as low as 10, and for a different quantity 
in another model too, as discussed in the next section. 
To rationalize this behavior one can resort to analytical 
approximations to the spectra. In particular, a semiclas- 
sical approach at Xc gives [19] 



Ek 



4/3 



which yields Tr 



fc^TV 1 / 3 and T C1 



(10) 



k- 1/3 N^\ k 
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III. THE DICKE MODEL 




FIG. 4. (Color online) Revival time as a function of \ f° r 
wave-packets centered around the ground state and TV = 1000. 
The inset is an enlargement of the delimited area and shows 
the divergence of Tr close to the critical point. Notice the 
logarithmic scale of the y axis. Time is given in dimensionless 
units. 



The Dicke model describes the interaction of a two- 
level atomic ensemble with a one-mode radiation field. 
It dates back to the 1950s [55] and to date several ex- 
perimental realizations of the model have been proposed 
|23) . The Dicke model exhibits a quantum phase transi- 
tion at zero temperature as embodied in the Hamiltonian 



H = w Q J z + wa)a + -^=(a f + a)(J+ + J_). (11) 

Here, J z and J are the usual angular momentum oper- 
ators for collective spin operators of length j = TV/2, 



and 



it 



are the bosonic operators of the field. The 



atomic level splitting is given by Wq, w is the field fre- 
quency, and A is the atom-field coupling. In the ther- 
modynamic limit, N,j — > oo, the system undergoes a 
second-order quantum phase transition at a critical cou- 
pling of A = A c = ^/wqw/2. At this point the sys- 
tem changes from a normal phase to the so-called super- 
radiant one in which both the field and the atomic en- 
semble acquire macroscopic occupations. Revivals and 
fractional revivals in the Jaynes-Cummings model, a sin- 
gle (j — 1/2) two-level atom interacting with a one mode 
of the quantized radiation field, were studied in the past 




FIG. 5. (Color online) Numerical spectra at criticality for in- 
creasing system-sizes TV as compared with reults of the semi- 
classical approximation. (Dimensionless units.) 



Proceeding along the same lines as in the previous 
section, that is, solving numerically for the spectra and 
populating the autostates around the ground state with 
Gaussian weights, a wave-packet is constructed and its 
time development monitored. The parity has been taken 
into account as a symmetry in this system. (Variational 
approximations have been proposed in [27J HE] to study 
the Dicke model in an analytical framework). We have 
verified that our results for the classical time on scale 
resonance, i.e. Wq = w = 1, follow those reported in |29j . 
A ~ (A c - A) 1 / 2 . As for T K , we find T R ~ |A C - Al" 1 (see 
Fig. [6]). Oddly enough, as in the vibron model, there 
is a value of A at which Xr diverges irrespective of the 
system-size. At this point, it may be in order to mention 
that it has been recently shown that the fidelity necessar- 
ily diverges at the critical point in the Dicke model [30] . 
and, interestingly, divergent behavior has been observed 
to occur at finite system-sizes [3"T] . 



IV. CONCLUSIONS 



being again the quantum number of the energy level the 
packet is centered around. However, this approximation 
seems to be valid only far from the ground state, as shown 
in Fig. [H]by plotting — Xc vs. k in double- logarithmic 
scale for several system sizes. Note that the spectra com- 
pare well with the power law k 136 only for sufficiently 
high k. This conveys the idea that the semiclassical for- 
mula might be asymptotically correct at very large values 
of both k and TV. 



We have extended the concept of quantum revivals 
to many-body systems and explored the implications of 
traversing a quantum phase transition by analyzing the 
time development of wave-packets centered around the 
fundamental state in two different models, namely, the 
U(3) vibron model and the Dicke model. Far from the 
quantum critical point, characteristic time scales such 
as the classical period and the revival and the super- 
revival times exhibit a smooth non-singular behavior, 
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FIG. 6. (Color online) Log-log plot of the revival times for 
the Dicke model with j = 10 in the vicinity of the critical 
point. The red circles and the blue squares correspond to, 
respectively, A < A c and A > A c . The data can approximated 
by a straight line of slope -1. (Atomic units.) 



packet with a Gaussian-distributed population basically 
a combination of states with almost equal energies. We 
have shown evidence that under these circumstances, the 
above-mentioned quantities diverge as power laws with 
well-defined critical exponents. Interestingly, the revival 
and super-revival times were found to diverge on ap- 
proaching the critical point even at finite system size. 
Should it be the case that all time scales in the Taylor 
expansion diverge, a statement we cannot confirm but 
that is suggested by the numerics and by the squeez- 
ing of the energy levels around the ground state as the 
transition point is approached, the time evolution of such 
wave packets at criticality would then be limited to phase 
changes, i.e., to rotations. Finally, we comment that 
quantum phase transitions can also influence the revival 
behavior of wave packets centered around excited states. 
This interesting effect, however, will be discussed else- 
where. 



ACKNOWLEDGMENTS 



and revival times are proportional to the system size. 
Upon approaching a quantum critical point, energy lev- 
els squeeze into the ground state, rendering the wave 



This work was supported by the Spanish Projects No. 
MICINN FIS2009-08451, No. FQM-02725 (Junta de An- 
dalucfa), No. 20F12.41 (CEI BioTic UGR), and No. 
MICINN FIS201 1-24149. 



[1] S. Sachdev, Quantum Phase Transitions (Cambridge 

University Press, New York, 1999). 
[2] R.W. Robinett, Phys. Rep. 392, 1 (2004). 
[3] E. Romera and F. de los Santos Phys. Rev. B 80, 165416 

(2009). 

[4] A. Lopez, Z.Z. Sun, and J. Schliemann, Phys. Rev. B 85, 
205428 (2012) 

[5] G. Rempe, H. Walther, and N. Klein, Phys. Rev. Lett. 

58, 353 (1987). 
[6] J. A. Yeazell, M. Mallalieu, and CR. Stroud Jr., Phys. 

Rev. Lett. 64, 2007 (1990). 
[7] T. Baumert, V. Engel, C. Rottgermann, W.T. Strunz, 

and G. Gerber, Chem. Phys. Lett. 191, 639 (1992). 
[8] M.J.J. Vrakking, D.M. Villeneuve, and A. Stolow, Phys. 

Rev. A 54, R37-R40 (1996). 
[9] A. Rudenko, Th. Ergler, B. Feuerstein, K. Zrost, CD. 
Schroter, R. Moshammer, and J. Ullrich, Chem. Phys. 
329, 193 (2006). 
[10] I.-Sh. Averbukh, M.J.J. Vrakking, D.M. Villeneuve, and 

A. Stolow, Phys. Rev. Lett. 77, 3518 (1996). 
[11] M. Mehring, K. Mueller, I.-Sh. Averbukh, W. Merkel, 
and W.P Schleich, Phys. Rev. Lett. 98, 120502 (2007). 
[12] E.A. Shapiro, M. Spanner, and M.Y. Ivanov, Phys. Rev. 

Lett. 91, 237901 (2003). 
[13] M. Spanner, E.A. Shapiro, and M.Y. Ivanov, Phys. Rev. 

Lett. 92, 093001 (2004). 
[14] K.F. Lee, D.M. Villeneuve, P.B. Corkum, and E.A. 
Shapiro, Phys. Rev. Lett. 93, 233601 (2004). 



[15] S. Will, T. Best, U. Schneider, L. Hackermuller, D.S. 
Luhmann, and I. Bloch, Nature (London) 465, 197 

(2010) . 

[16] E. Tiesinga and PR. Johnson, Phys. Rev. A 83, 063609 

(2011) . 

[17] F. Iachello and S. Oss, J. Chem. Phys. 104, 6956 (1996). 
[18] F. Perez-Bernal and F. Iachello, Phys. Rev. A 77, 032115 
(2008). 

[19] M.A. Caprio, P. Cejnar, and F. Iachello, Ann. Phys. 

(NY.) 323, 1106-1135 (2008). 
[20] In principle, a should be chosen such that the energy 

variances for different N are comparable, but we have 

checked that this is not a relevant issue. 
[21] P. Perez-Fernandez, J.M. Arias, J.E. Garcfa-Ramos, and 

F. Perez-Bernal, Phys. Rev. A 83, 062125 (2011). 
[22] R.H. Dicke, Phys. Rev. 93, 99 (1954). 
[23] B.M. Garraway, Philos. Trans. R. Soc. A 369, 1137-1155 

(2011). 

[24] C. Emary and T. Brandes, Phys. Rev. E 67, 066203 
(2003). 

[25] J. Gea-Banacloche, Phys. Rev. Lett. 65, 3385 (1990). 

[26] I.-Sh. Averbukh, Phys. Rev. A 46, R2205 (1992). 

[27] O. Castahos, E. Nahmad-Achar, R. Lopez-Pena, and 

J.G. Hirsch, Phys. Rev. A 83, 051601 (2011). 
[28] O. Castahos, E. Nahmad-Achar, R. Lopez-Pena, and 

J.G. Hirsch, Phys. Rev. A 84, 013819 (2011). 
[29] J. Vidal and S. Dusuel, Europhys. Lett. 74, 817 (2006). 
[30] O. Castahos, E. Nahmad-Achar, R. Lopez-Pena, and 

J.G. Hirsch, Phys. Rev. A 86, 023814 (2012). 



6 



[31] W.G. Wang, P. Qin, Q. Wang, G. Benenti, and G. Casati 
Phys. Rev. E 86, 021124 (2012). 



